宽 视 场 成 像 网 格 化 算法 中 w-plane 最 优 经 验 值 研究 
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摘要 : 射电 干涉 阵列 宽 视 场 成 像 网 格 化 处 理 过 程 中 必须 考虑 w 项 的 影响 。w-projection f w- 
stacking 是 两 个 重要 的 宽 视 场 成 像 网 格 化 处 理 算 法 ， w-plane 参数 是 算法 中 影响 计算 速度 和 成 图 质量 的 
一 个 重要 因素 。 本 文 研究 了 w-projection 和 w-stacking 两 种 网 格 化 算法 ， 利 用 SKA1 低频 阵 台 站 数据 和 
ASKAP 软件 包 进 行 模拟 观测 ， 对 两 种 算法 在 不 同 w-plane 参数 取 值 情况 下 的 成 图 速度 和 成 图 质量 进行 了 定 
量 分 析 对 比 。 nea eid w-plane 是 性 能 改善 的 重要 参数 ， 针 对 woprojection 算法 来 说 ，w-plane 
i 经 验 值 要 大 才能 较 好 的 得 到 成 像 效 果 。w-stacking 算法 虽然 有 很 大 的 速度 优势 ， 但 
算法 实现 中 w-plane ee 论文 给 出 了 推荐 的 w-plane 取 值 。 本 文 工 作 是 大 视 场 成 像 算法 的 
基础 研究 工作 ， 对 未 来 SKA 科学 数据 处 理 中 的 管线 设计 有 较 好 的 参考 价值 。 
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中 图 分 类 号 : P161 文献 标识 码 : A ”文章 编号 : 1672-7673(2018)xx-xxxx-xx 


平方 公里 阵 SKA (square kilometer array) 将 是 国际 上 建造 的 最 大 综合 孔径 射电 望远镜 中。SKA 由 数 
于 个 反射 面 天 线 和 多 达 一 百 万 个 低频 天 线 组 成 ， 总 接收 面积 达到 一 平方 公里 。 

SKA 科学 数据 处 理 (Science Data Processor, SDP) 是 SKA 设计 中 的 一 个 重要 组 成 部 分 。SDP 专注 于 设 
计 硬 件 平台 、 软 件 和 算法 来 处 理 射 电 望 远 镜 的 观测 数据 户 。SDP 的 设计 面临 的 挑战 是 多 方面 的 。SDP 需 
要 处 理 来 自 望远镜 的 海量 观测 数据 量 。SKA1I1( 约 占 10% 的 总 体 规 模 ) 数 据 量 将 达到 300PB/ 年 。 为 了 满足 
SKA 科学 研究 开展 的 需要 ，SDP 在 数据 处 理 管 线 设 计 中 ， 需 要 充分 地 考虑 数据 处 理 的 性 能 与 质量 。 由 于 
SDP 处 理 需 要 的 计算 力 高 出 现 有 最 大 的 天 文 数据 处 理 系统 两 个 数量 级 ，SDP 的 关键 算法 的 运行 速率 要 达 
到 百 万 兆 级 的 运算 水 平 。SKA2 的 计算 力 需 求 估计 需要 达到 4 Exaflops 的 水 平 。 因 此 ， 开 展 对 SDP 数据 
处 理 中 关键 算法 的 研究 ， 提 高 数据 处 理 的 性 能 具有 显著 的 意义 。 

网 格 化 《Gridding) 算法 是 综合 孔径 望远镜 数据 处 理 中 的 一 个 重要 算法 ， 对 计算 的 消耗 占据 了 很 大 比 
重 。 对 网 格 化 算法 进行 性 能 优化 ， 对 SDP 有 显著 的 价值 和 意义 。 低 频 射 电 望 远 镜 大 视 场 成 像 受到 很 多 条 
件 的 制约 ， 其 中 最 重要 的 一 个 影响 因素 是 非 共 面 基线 效应 。 在 大 视 场 与 非 共 面 基线 效应 下 ， 对 可 见 度数 
据 直 接 进 行 傅立叶 变换 会 导致 图 像 严重 蚁 变 ， 因此 ， 在 大 视 场 成 图 时 ， 解 决 w 项 的 问题 至 关 重 要 。 目 
前 ， 应 对 大 视 场 与 非 共 面 基线 效应 的 算法 有 : facetingBl]、 三 维 傅立叶 变换 内 、w-projection5]、w-stackingtg 
和 warped snapshots[41 等 算法 以 及 其 他 复合 算法 。 在 这 些 方 法 中 ，w-projection 和 w-stacking 是 目前 应 用 相 
对 较 广 的 网 格 化 算法 。 其 中 w-projection 算法 在 牺牲 较 大 内 存 的 情况 下 ， 具 有 优越 的 计算 速度 和 误差 控 
制 。 该 方法 不 仅 可 以 使 用 传统 的 洁 化 (CLEAN) 方 法 去 卷 积 ， 还 可 以 选择 多 尺度 洁 化 方法 ， 在 高 动态 范围 
下 ， 有 较 大 的 优势 , 相位 中 心 基本 不 失真 。w-stacking 算法 在 图 像 空间 以 消耗 更 多 的 内 存 为 代价 处 理 傅 里 
叶 空 间 中 的 w 项 ， 计 算 速 度 优 于 w-projection 快门。 
在 w-projection 和 w-stacking 算法 的 具体 实现 过 程 中 ，w-plane 是 一 个 重要 参数 。 在 w 的 最 大 值 和 最 
小 值 之 间 分 出 多 个 w-plane， 对 不 同 w-plane 的 可 见 度数 据 使 用 不 同 的 卷 积 核 进行 卷 积 运算 。 因 此 ，w- 
plane 的 数量 直接 影响 到 算法 的 计算 速度 、 内 存 开销 以 及 成 图 质量 。 本 文 主要 对 w-projection 和 w-stacking 
算法 在 实现 与 实际 运行 中 的 w-plane 进行 实验 研究 ， 利 用 SKA1 低频 阵 台 站 数据 和 ASKAP 软件 包 进 行 模 
拟 观测 ， 以 期 掌握 w-plane 对 整体 计算 性 能 的 影响 ， 为 后 期 的 SKA-SDP 建设 以 及 国内 科学 数据 中 心 建设 
中 的 管线 系统 设计 提供 参考 。 
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下 ow 项 的 值 将 远 远 大 于 1， 因 此 直接 对 可 见 度数 
题 ， 天 文学 家 提出 了 


w-plane 


低频 射电 望远镜 大 视 场 成 像 最 重要 的 


RING 
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影响 因 


进行 二 维 健生 
法 用 来 应 对 非 共 面 基 线 效应 造成 的 图 


素 就 是 非 共 面 基线 效应 ， 在 非 共 面 基 线 效 应 的 影响 
有 叶 变 换 的 条 件 已 经 不 满足 ， 针 对 这 一 问 
像 畸变 。w-projection 算法 通过 将 不 同 


w 值 的 可 见 度 数据 投影 到 w=0 平面 ， 从 而 消除 w 项 的 影响 ， 有 效 提高 成 图 质量 。w-stacking 算法 能 够 在 


ica 


F 面 矫正 非 共 面 
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起 的 相位 人 


1.1 w-projection 算法 中 的 w-plane 


度数 据 ， 在 计算 过 程 中 应 该 采 
大 和 最 小 值 之 | 
对 每 个 可 见 度数 据 与 


在 射电 综合 孔径 成 像 处 型 


Viu, v, w) = [. 全 


w-projection 算法 将 不 同 w 值 的 可 见 度数 据 投 影 到 w=0 平面 ， 


Viu,v, w = 0) = [. T 


I(l, m) 


y1- P -m° 


ij 移 从 而 有 效 还 原 天 空 亮 度 分 布 。 


eilutvm tr wy -m° -dm 


(4) 


I(l, m) 


VI- -m 


经 过 简单 的 指数 运算 ， 式 (1) 变化 为 式 (3) 和 (4) : 


vvn) = EE 


公式 (3) 、 


上 述 过 程 显 示 ， 在 w-projection HIE, ARREST 
不 同 的 卷 积 核 。 在 实 
本 平均 分 成 多 个 w-plane， 然 后 对 每 个 w-pane 按照 公式 (6) 计生 
其 距离 最 近 的 w-plane 的 卷 积 核 进行 卷 积 ， 从 而 实现 到 w=0 平面 的 投影 。 
量 和 内 存 消 耗 大 ，w-plane 的 值 
w-plane 的 数量 为 零 ， 则 直接 退化 为 二 维 傅立叶 变换 。Cornwell 


显而易见 ，w-plane 的 值 
量 和 内 存 消耗 小 ， 但 是 投影 效果 欠 佳 。 当 


(4) 可 以 改写 为 公式 (5) ， 其 
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法 的 核心 。 到 


公式 C1) 简化 为 (2): 


e= +m) \dldm 
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卷 积 核 通过 公式 (6) 近似 计算 得 到 : 


际 计算 过 程 中 ， 


R, Biri 


等 在 前 期 研究 中 取 最 大 的 w 值 的 平方 根 作为 经 验 值 品 。 


2 w-stacking 算法 中 w-plane 


w-stacking 算法 的 提出 也 是 为 了 消除 非 共 面 基线 效应 对 大 视 场 成 像 的 影响 ， 但 是 不 同 于 w-projection 算 
法 ，w-stacking 并 没有 消除 w 项 ， 而 是 对 其 进行 了 修正 。 


根据 公式 (1) 将 uv 空间 转换 到 Im 空间 得 到 式 〈7) : 
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过 程 中 ， 可 见 度 数据 与 天 空 亮度 之 间 的 关系 如 公式 〈1) 所 示 : 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


E 想 情况 下 ， 位 于 不 同 w-plane 的 可 见 
般 计 算 w 的 最 大 和 最 小 值 ， 在 最 
保存 下 来 ， 最 后 


越 小 ， 计 算 


(7) 


对 公式 (7) 的 两 边 在 w 范围 内 进行 积分 得 到 式 C8) : 

I, m) (Wa E Wain) ia eons Bay? D V (u, v, wel) dudvdw (8) 

VI-P -mw Pain -o J-e 
其 中 wy 和 wi, 代表 w 的 最 大 值 和 最 小 值 。 式 (8) 的 等 号 右 侧 包 括 两 部 分 ， 一 是 二 维 傅立叶 变换 ， 二 
是 对 w 的 复数 积分 。 对 式 (8) 的 等 号 右 侧 离散 化 得 到 式 CO) : 
I(l, m) (Wax = Win) _ > gO- -Do Wu, v, w,)} (9) 
Ji = P = m’ a=Wnin 

see FW Cu, v, w, 间 表示 将 可 见 度数 据 在 特定 的 离散 w-plane 上 进行 网 格 化 。 
w-stacking 算法 的 实现 步骤 如 下 ; 

(1) 根据 可 见 度 数据 w 的 范围 设 定 w-plane 的 细 分 粒度 ; 

(2) 依据 每 个 可 见 度 数据 的 w 值 ， 将 可 见 度 数据 网 格 化 到 某 一 w-plane; 

G) 对 每 个 w-plane 的 可 见 度数 据 进行 二 维 傅立叶 变换 得 到 相应 图 像 ; 

(4) 根据 w 值 对 各 个 w-plane 的 图 像 进行 相位 矫正 ; 

(5) 加 权 有 登 加 所 有 图 像 ， 得 到 最 终 图像 。 
在 w-stacking 算法 实现 过 程 中 ，w-plane 分 的 过 粗 会 导致 产生 伪 影 ， 分 的 过 细 又 会 增 大 计算 量 和 内 存 消 


耗 


同时 ，w-stacking 算法 将 不 同 w-plane 上 的 可 见 度数 据 看 成 是 相互 独立 ， 各 个 面 的 数据 可 以 并 行进 行 
二 维 傅立叶 变换 和 相位 改正 ， 因 此 ，w-stacking 算法 在 并 行 实现 方面 具有 优势 。 


3 观测 模拟 仿真 实验 


3.1 实验 环境 


为 了 验证 w-plane 的 影响 ， 以 期 获得 较为 理想 的 w-plane 取 值 ， 在 确保 成 像 质量 的 情况 下 降低 成 像 处 
理 时 间 ， 本 文 进 行 了 一 系列 实验 。 数 据 处 理 软件 采用 了 澳大利亚 SKA 探 路 者 (Australian SKA Pathfinder, 
ASKAP) 的 数据 处 理 软件 系统 四 。ASKAP 通过 名 为 ASKAPsoft 的 软件 系统 进行 望远镜 观测 控制 和 数据 处 理 
Ol, ASKAPsoft 框架 由 三 部 分 组 成 : 望远镜 观测 系统 、 处 理 中 心 和 科学 数据 存档 。ASKAP 通过 bash shell 
脚本 组 合 ASKAPsoft 中 的 程序 形成 数据 处 理 管线 ， 主 要 的 数据 处 理 管线 有 数据 摄取 管线 、 定 标 管线 和 成 
图 管线 00。 

本 实验 在 8 台 汉 柏 服务 器 上 完成 ， 有 具体 硬件 配置 如 表 1 所 示 。 

表 1 硬件 设备 参数 


Table 1 The parameters of device 


硬件 设备 设备 参数 

服务 器 型 号 C640 G3 (8) 

CPU Intel Xeon E5-2620 v3 

内 存 128GB 

网 络 Infiniband + 10Gb 以 太 网 


3. 2 模拟 观测 与 成 像 处 理 


为 了 分 析 w-plane 对 w-projection 和 w-stacking 中 的 影响 ， 我 们 采用 了 模拟 观测 、 成 像 、 对 后 续 图 像 进 
行 分 析 的 方法 来 验证 不 同 的 w-plane 造成 的 影响 。 

1、 模拟 观测 

实验 采用 SKA1-Low 望远镜 的 512 个 天 线 台 站 坐标 数据 ， 天 线 布局 如 图 1 所 示 。 天 空 模型 由 73 个 点 
源 构成 ， 整 体形 状 为 “ 米 ” 字 ， 如 图 2 (a) 所 示 。 实 验 采 用 16 个 计算 节点 ， 每 个 计算 节点 处 理 一 个 通道 的 数 
据 ， 通 道 带 宽 1IMHZ， 因 此 观测 频率 履 盖 范围 16MHz， 积 分 时 间 间 隔 150s， 模 拟 观 测 时 长 6 小 时 。 


SKA1-Low 天 线 频率 履 盖 范围 50-350MHz， 本 实验 中 心 参考 频率 为 100MHz， 设 置 16 个 带宽 为 1IMHz 的 
通道 ， 因 此 观测 频率 范围 为 92-108MHz。 
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图 1 天 线 分 布 
Fig. 1 Antennas layout 
使 用 模拟 观测 程序 对 天 空 模型 进行 16 个 通道 的 模拟 观测 ， 得 到 18837504 行 5.5GB 的 模拟 可 见 度数 


据 。 

2、 成 像 处 理 
使 用 w-projection 和 w-stacking 算法 ， 对 模拟 观测 获得 的 可 见 度数 据 16 个 通道 进行 连续 谱 成 图 。 
ASKAP 成 图 管线 在 8 个 计算 节点 上 并 行 运算 ， 每 个 计算 节点 处 理 两 个 通道 。 实 验 过 程 中 ， 不 断 改变 w- 
plane 的 数量 ， 得 到 不 同 w-plane 设置 下 的 图 像 ， 并 记录 每 次 成 图 时 间 供 后 续 分 析 使 用 。 图 4 (a) 显示 了 
随 着 w-plane 数量 的 增长 ， 两 种 算法 处 理 时 间 的 变化 。 

图 2 (b) 是 二 维 傅立叶 处 理 得 到 的 图 像 ， 可 以 看 到 w 项 影响 导致 的 拖 尾 效应 (Smearing) 非 常 明显 。 图 3 
是 w-plane 数量 设置 为 171 的 情况 下 ，w-projection 和 w-stacking 算法 处 理 得 到 的 图 像 ， 拖 尾 效应 得 到 抑 
制 ， 但 边缘 的 模拟 星 并 没有 被 很 好 地 还 原 。 
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图 2 (a) 天 空 模型 图 ; (b) 二 维 傅 里 叶 变 换 成 像 结果 


Fig. 2 (a) Sky model image; (b) The result of two-dimensional Fourier transform imaging 
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3、 图 像 质 量 
ASH 


Duchamp 星 源 搜索 算法 在 ASKAP 软件 
时 已 经 预 设 了 总 的 星 数 为 73， 同 时 所 有 星 的 位 置 已 知 。 
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评价 方法 
F 价 不 同 w-plane 数量 


=r 


下 的 两 和 


EA 


还原 的 星 数 量 


SELAVY 搜索 出 点 源 数量 的 变化 。 
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图 3 (a) w-projection 算法 的 成 图 结果 ; (b) w-stacking 算法 的 成 图 结果 
Fig. 3 (a) The result of the w-projection imaging; (b) The result of the w-stack imaging 
' 算 法 的 成 像 质量 ， 研 究 采 用 了 SELAVY 软件 包 ， 该 软件 包 是 
! 的 实现 模块 ， 用 来 搜索 观测 结果 中 的 星 数 量 。 由 于 在 仿真 观测 
因此 ， 通 过 SELAVY 搜索 出 的 星 与 已 知 天 空 模型 
量 来 快速 判断 成 像 的 质量 。 图 4 (b) 显示 了 随 着 w-plane 数量 的 增 
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Fig. 4 (a) The relation between the running time and w-plane; (b) The relation between the number of searched sources and w-plane 


3、 结 果 分 析 


本 实验 最 大 基线 长 度 81602 米 ， 最 大 w 约 为 28000 个 波长 。 根 据 
plane 值 一 般 设置 为 最 大 w 的 
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4 (a) 可 以 看 到 w-projection 的 运行 时 间 随 着 w-plane 值 的 增 


小 于 Vw 的 什 


后 其 增长 速度 很 
像 成 像 时 太 多 的 w-plane 需要 


w- faeries 算法 也 就 失去 了 速度 优势 。 


SELAVY 软件 包 对 两 
plane 较 小 的 情况 下 w-stacking 算法 成 
图 质量 近似 相同 。 
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期 研究 的 结果 ，w- 


的 情况 下 其 计 


得 到 星 的 数 
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生 增 长 。 而 w-stacking 算 
间 明 显 小 于 
F 1.5 Vw 之 后 超过 $ w-projection 算法 的 运算 
进行 相位 的 改正 从 而 导致 w-stacking 算法 的 计算 力 消 耗 


算法 在 不 同 的 w-plane 下 成 图 结果 进行 星 像 统计 。 从 图 4 (b) 可 以 看 到 在 w- 


w-projection 算法 ; 


E% F w-projection 的 结果 ， 但 是 随 着 w-plane 的 增加 


w-plane 的 取 值 不 同 于 以 往 的 经 验 值 Vw ， 实 验 可 以 看 出 : 在 wplane 大 于 1.5 Vw 之 后 成 图 质量 更 
好 ，SELAVY 程序 包 可 以 得 到 所 有 的 73 个 星 源 的 数据 ， 因 此 为 获得 高 分 辩 率 的 天 空 图 w-plane 的 取 值 应 
大 于 1.5 Vw 。 虽 然 两 者 的 成 图 质量 相同 但 是 w-stacking 算法 在 w-plane 的 值 小 于 Vw 的 值 的 情况 下 相 较 
于 w-projection 具有 很 大 的 速度 优势 ， 但 当 w-plane 的 值 大 于 Vw 后 ，w-stacking 的 速度 优势 将 不 存在 。 
通过 分 析 这 两 种 算法 的 成 图 质量 和 运行 时 间 趋 势 ， 可 以 看 到 w-stacking 算法 更 适合 w-plane 值 小 于 
Vw 的 情况 ， 为 弥补 其 weplane 较 小 引起 的 成 像 准 确 度 的 问题 ， 今 后 可 以 考虑 通过 结合 多 波束 观测 进行 成 
图 ， 再 通过 SKA 的 LINMOS(Linear Mosaic Applicaton 软 件 包 实现 多 波束 成 图 后 的 拼接 ， 该 软件 包 实 现 一 


系列 图 片 的 线性 拼接 ， 从 而 提高 成 像 质量 。w-projection 算法 则 更 适合 通过 提高 w-plane 的 值 实现 更 为 准 
确 的 成 图 o 


结论 


综 上 所 述 ， 针 对 w-projection 算法 和 w-stacking 算法 在 管线 中 的 应 用 ， 有 如 下 结论 
1、 合 理 的 w-plane 的 取 值 对 于 成 像 有 关键 的 作用 。 实 验 表 明 ，w-plane 的 取 值 应 是 经 验 取 值 Vw 的 
1.5 倍 ， 此 时 w-projection 算法 和 w-stacking 算法 都 可 以 得 到 很 好 的 成 图 质量 ， 这 与 Cornwell 估计 的 有 
i es 
w-plane 的 取 值 对 于 结果 的 正确 性 有 紧密 的 关系 。 事 实 上 只 有 w-plane 取 值 较 大 时 才 有 可 能 完整 地 
Peete SE A. ARR TEACHER EE, 2 CE HE Wh SE BE OR A EN Tad fg Ay TE AT TE. “4 w-plane 的 值 小 
于 Vw 的 时 候 ， w-stacking 算法 比 w-projection 算法 有 很 大 的 计算 时 间 优 势 。 对 于 成 图 分 辩 率 要 求 高 且 w- 
plane 的 值 大 于 Vw 的 时 候 ，w-stacking 算法 将 失去 速度 优势 。 
本 文 的 工作 表明 SKA-SDP 数据 处 理 管线 应 根据 不 同 的 需求 选择 不 同 的 ai 算法 ， 后 续 工 作 可 以 
在 运算 速度 提高 ， 提 高 算法 的 并 行 性 为 研究 重点 ， 以 期 获得 更 好 的 处 理 能 
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A Study on Optimal Experiential W-plane for Wide-Field 
Gridding Algorithm 


Yu Xiaoyul Deng Hui!?, Mei Ying!, Wei Shoulin!, Shi Congming!, Wang Feng!” 
(1. Key Laboratory of Computer Technology Application; 2. Center for Astrophysics, Guangzhou University) 


Abstract: Astronomy data processing platform is facing unprecedented big data challenge. SK A(square 
kilometer array) as the biggest radio synthesis telescope, its generated data rate of averaged visibilities will be 
hundred gigabit per second. SKA-SDP handles big data and extremely large computing cost by data processing 
pipeline. Gridding algorithm will take a large proportion of compute load. How to save power cost and improve 
the computational efficiency, on the prerequisite of guaranteed imaging quality, need to be well-considered. W- 


projection and w-stacking are two gridding algorithms possibly be used in SKA-SDP; The present paper study on 
the effect of w-plane in those algorithms. Based on SKA continuum imaging pipeline and simulation observation, 
we analyze the effect of w-plane to the imaging results. Our experimental results show that w-plane is a crucial 
parameter for performance promotion, and w-plane should be bigger than experimental value in w-projection 
algorithm. Compare with w-projection algorithm, w-stacking have advantage in run time, yet w-plane have great 
effect on w-stacking algorithm. The present paper is a basic work for study on large field imaging, and a new w- 
plane for w-projection and w-stacking algorithm is proposed. Our study can serve as a reference for designs of 
SKA data processing pipeline. 

Key words: SKA-SDP; Wide-field imaging; Gridding algorithm 


